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INTRODUCTION 

Due to the wide variety of uses of sensitivity derivatives, the development of 
efficient computational procedures for the calculation of these derivatives has at- 
tracted considerable attention in recent years. The calculation of sensitivity 
derivatives forms the backbone of many optimization procedures and is the major 
contributor to the cost and time of optimization of large systems. In addition, 
sensitivity derivatives have several other applications in structural mechanics in- 
cluding approximate analysis (and reanalysis) techniques, analytical model improve- 
ment, and assessment of design trends. A review of the state of the art in 
sensitivity calculations is contained in a survey paper (Ref. 1), a monograph 
(Ref. 2), as well as in some papers in these proceedings. Despite all the recent 
advances made, the calculation of sensitivity derivatives for large structural 
systems (with large number of degrees of freedom and design variables) , is quite 
expensive even on present-day large computers. 

The present study focuses on the development of efficient techniques for cal- 
culating sensitivity derivatives. Specifically, the objective and scope of the 
present paper are listed in Fig. 1. The objective is to present a computational 
procedure for calculating sensitivity derivatives as part of performing structural 
reanalysis for large-scale problems. The scope of the paper is limited to framed 
type structures. Both linear static analysis and f ree-vibration eigenvalue problems 
are considered. 


Objective 

• To present a computational procedure for calculating 
sensitivity derivatives as part of performing structural 
reanalysis for large-scale problems 


Scope 


• Fra me -type structures 

• Linear static analysis 

• Eigenvalue problems 

Figure 1 
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BASIC IDEA AND KEY ELEMENTS OF THE PROPOSED PROCEDURE 


The basic idea and the three key elements of the proposed procedure are listed 
in Fig. 2. The basic idea is to generate the solution of the modified structure 
using large perturbations from that of the original structure. The three key ele- 
ments are: a) lumping of the large number of design variables into one (or a small 
number of) tracing parameters; b) application of operator splitting/reduction 
technique; and c) for very large problems, use of single-level or multilevel sub- 
structuring. Only the first two key elements are discussed in this paper. The 
application of operator splitting/reduction technique proved to be effective in 
reducing the computational effort in a number of structural mechanics problems (see, 
for example. Refs. 3 to 6). 


Basic idea 

• Solution for modified structure is obtained using large 
perturbations from that of original structure 

Key elements 

• Lumping of design variables into tracing parameter(s) 

• Application of operator splitting/ reduction technique 

• Use of multilevel substructuring (for very large problems) 


Figure 2 
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APPLICATION TO LINEAR STATIC ANALYSIS 


In Fig. 3 the application of the proposed procedure to linear static analysis 
is outlined. The governing finite element equations of the original and modified 
structures are shown. The global stiffness matrices, load vectors, and responses 
of the original and modified structures are designated by [k] q , [k] ; {P}; and 

{X} , {X}, respectively. The original and modified structure characteristics cor- 
^ o * 

respond to the values of d^ and d^ of the design variables, respectively. 

The operator splitting technique is now applied, and the equations of the 
modified structure are expressed in terms of the original structure equations plus 
correction terms. A tracing parameter X is introduced and is attached to the cor- 
rection terms. The tracing parameter is dimensionless and identifies alt the design 
modifications. The original structure equations correspond to A=0, and the modified 
structure equations correspond to A=l. 


Operator splitting 
Original structure 

MoiX'o = iP, o : Wo=[ K,d Mo 

Modified structure 

# 

[k] 1X1= {Pt; [K] = [Kid,)] 

'[M„ 

X = 0 — Original structure 
X = 1 — Modified structure 

Figure 3 
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REDUCTION METHOD FOR STATIC PROBLEMS 


The response of the modified structure, {X}, is now expressed as a linear com- 
bination of a few preselected global approximation vectors (or modes). This is 
expressed by the transformation shown in Fig. 4. The columns of the matrix [r] are 
the global approximation vectors, and the elements of the vector are the ampli- 
tudes of the approximation vectors which are, as yet, unknowns. Note that the 
number of global approximation vectors, r, is considerably smaller than the total 
number of degrees of freedom, n. 

A Rayleigh-Ritz technique is now used to approximate the governing equations of 
the modified structure by a much smaller system of equations in the unknowns {^}. 


Basis reduction 



; r < < n 


where {if/} = amplitudes of global approximation vectors 


Reduced system of equations 

Rayleigh -Ritz technique used to approximate the equations of 
the modified structure 

[M 0 + X (M - 

where [k] q = |r]*[K] o [r| 

[k] = [r|*[K|Ir| 

n = iifn 

)Pl = (rlV) 


Figure 4 
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SELECTION AND EVALUATION OF THE GLOBAL APPROXIMATION 


The effectiveness of the proposed procedure depends, to a great extent, on the 
proper choice of the global approximation vectors. In the present study the global 
approximation vectors are selected to be the response of the original structure, 

{X} q , and its various-order derivatives with respect to the parameter X. The re- 
cursion relations for evaluating the approximation vectors are obtained by successive 
differentiation of the original finite element equations. Note that the matrix on 
the left hand sides of these equations, [k] , is the same (see Fig. 5). 


fxt = [r] w 

where [r] = [W 0 l ^( 0 . .] 

0 A 

Mo ,Xl o = {P, o 

M 0 {f| o = -(M-[ k ] 0 ) M o + |P( - 1P| o 

- 2 (m 


Same left hand side. 

Figure 5 
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COMPUTATIONAL PROCEDURE 


The computational procedure for generating the solution of the modified struc- 
ture and the sensitivity derivatives is outlined in Fig. 6. 

The first step is to generate the global approximation vectors at A=0 through 
decomposition of the stiffness matrix of the original structure. The derivatives 
with respect to X provide information about the sensitivity of the response to dtt 
the design modifications. Because X is dimensionless, the derivatives with respect 
to X have the same dimensions as the original response quantities, and consequently 
an assessment of the effect of design modifications on the response can be easily 
made. 


The second step is to generate the reduced equations and solve them for the 
amplitudes of the global approximation vectors. 


• Evaluate global approximation vectors at A = 0 

• Derivatives with respect to A represent sensitivity of 
the response to design modifications 


• Generate reduced equations 


• Solve reduced equations and find amplitudes of 
global approximation vectors 


Figure 6 
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RELATIONSHIP BETWEEN THE PRECONDITIONED CONJUGATE GRADIENT (PCG) 

TECHNIQUE AND THE PROPOSED COMPUTATIONAL PROCEDURE 

If the proposed computational procedure is contrasted with the preconditioned 
conjugate gradient (PCG) technique in which the preconditioning matrix is selected 
to be the global stiffness matrix of the original structure, [k] , the relationships 
shown in Fig. 7 can be identified. These relations express the preconditioned re- 
siduals {y} Q , {y}^, ... of the PCG technique in terms of the global approximation 

2 

vectors of the foregoing technique, {— } , { — -} , 

8A 0 ax 2 o 


Equivalence 

[Wo + KM - M 0 )] ,xl= « p > 0 + *(< p > - < p > 0 ) 


[k] = preconditioning matrix 


PCG 

Proposed procedure 

Wo — 

Preconditioned 

residuals 

— {X} 

0 

\ 

Mo — 

/ vl 

(dX 

(dX 


i y h 

Mi — 

(1 V iaxj 

— ¥ c. 

i=i 1 

. 2 LM 

0 0 A 0 

\*k\ 

0 


Figure 7 
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IMPLICATIONS OF THE SIMILARITIES BETWEEN THE PROPOSED 
PROCEDURE AND THE PCG TECHNIQUE 

The implications of the similarities between the proposed computational pro- 
cedure and the PCG technique are listed in Fig. 8. 

For the PCG technique, the similarities can be exploited to provide a rational 
approach for selecting the preconditioning matrix (as the global stiffness matrix 
of the original structure), and a physical meaning for the preconditioned residual 
vectors (in terms of sensitivity derivatives) . 

For the proposed procedure, some of the work done on parallelizing the PCG on 
multiprocessor computers can be exploited. 


• Rational choice for preconditioning matrix 

• Physical meaning for preconditioned residuals 
(in terms of sensitivity derivatives) 


Proposed procedure 

• Exploiting work done on parallelizing PCG on 
multiprocessor computers 


Figure 8 
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APPLICATION TO EIGENVALUE PROBLEMS 


The application of the proposed computational procedure to free vibration 
(eigenvalue) problems is outlined in Fig. 9. The governing equations of the origin- 
al and modified structures are shown. Again, the operator splitting technique is 
applied, and the stiffness and mass matrices of the modified structure are written 
as the sum of the corresponding matrices of the original structure plus correction 
terms. The correction terms are identified by the tracing parameter X. When X=0, 
the original structure equations are recovered, and when X=1 the modified structure 
equations are obtained. 


Operator splitting 

Original structure 

[Wo - fi Mo]' X 'o= 0 
Modified structure 

f[K] - fi[lVl]]|X(= 0 

or, 

[(Mo - fi M 0 ) + *(M a - q[m] 8 )]m= 0 

A = 0 ► original structure 

X = 1 ► modified structure 

where [k] a = [k] - [k] q 


Figure 9 
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REDUCTION METHOD FOR EIGENVALUE PROBLEMS 


The application of the reduction method to the free vibration problem is outlined 
in Fig. 10. As in the static analysis, the eigenvectors of the modified structure, 
{X}, are approximated by a linear combination of a few global approximation vectors. 
This is accomplished by the transformation shown. An efficient choice of the ap- 
proximation vectors was found to be a few eigenvectors for the original structure 
(corresponding to A=0) and their derivatives with respect to A, evaluated at A=0. 

Then, the Rayleigh-Ritz technique is used to approximate the original large 
eigenvalue problem by the reduced one shown in Fig. 10. The solution of the reduced 
eigenvalue problem gives the amplitudes of the global approximation vectors. 


Basis reduction 


• Eigenvectors of 
IX} 

n 


modified structure, {X|, are approximated by ; 
1= Mn.r^'r.1 ;r<< " 


where 





J A = 0 


Reduced system of equations 


• Rayleigh-Ritz technique is used to approximate the original 
eigenvalue problem by a reduced one 


where 




+ *(Ma ' Q Ma) 


101=0 


[i<l 0 = I r I* l K l o l r l 
|i<l a = l r l* l K l a |r) 

No = l r l'l M l 0 |r| 

M a = l r l t M a l r l 


Figure 10 
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EVALUATION OF GLOBAL APPROXIMATION VECTORS 


The recursion equations used in generating the eigenvectors at A=0 and their 
first derivatives with respect to A are shown in Fig. 11. Note that the left hand 
sides of all these equations are the same. The expression of the first derivatives 
of the eigenvalues with respect to A appearing on the right hand sides of the equ- 
ations are given in Fig. 11. 

Since the matrix on the left hand side of the recursion equations used in 
evaluating these derivatives is singular, the solution of each set of equations can 
be expressed as the sum of a homogeneous solution (multiple of the eigenvector) and 
a particular solution, {Q}^. The equations used in evaluating the particular solu- 
tion {Q} are given in Fig. 11. The details of this procedure are given in Ref. 5. 


Recursion formulas 
Jk] 0 - G [m] J (XI = 0 

Mo' Q Mo]i^| = df [Wo + X Ma] )x, -[[ K ]a • fi [ M ] a ] ,X| 
tf = < x »*[M a -«N a ]t x l/fXl t [[M] o+ x[M] a ]|X| 


{o| T ,Q >1 + C 1 w 

Where {Qij = particular solution. 


Figure 11 
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COMPUTATIONAL PROCEDURE FOR EIGENVALUE PROBLEMS 


The procedure for extracting the eigenvectors of the modified structure and for 
generating the sensitivity of the eigenvectors to design modifications is outlined 
in Fig. 12. 

First: A few eigenvectors of the original structure (corresponding to A=0) are 

generated . 

Second: The derivatives of the eigenvectors with respect to A are generated at 

A^O. In the process, derivatives of the eigenvalues are also computed. These de- 
rivatives provide sensitivity information regarding the effect of all the design 
modifications on the eigenvectors and eigenvalues. The reduced equations are 
generated . 

Third: The reduced eigenvalue problem is solved at A=l. 


• Generate eigenvectors for original structure (A = 0) 


• Generate global approximation vectors (derivatives of 
eigenvectors w. r.t. A) and reduced equations. In the 
process, derivatives of Qw.r.t. A are computed 


• Solve reduced eigenvalue problem at A = 1 


Figure 12 
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CANTILEVERED LATTICE TRUSS 


To assess the effectiveness of the proposed computational procedure, a number 
of problems were solved by this procedure. Comparison was made with the direct solu- 
tion of the structure. Herein a typical problem of a five-bay cantilevered lattice 
truss is considered (see Fig. 13). In the original structure all the longerons had 
the same cross section, and all the battens and diagonals had the same cross section. 
The design variables consisted of the cross sectional areas, moments of inertia and 
torsional constants. The characteristics of the original and modified structures 
are given in Fig. 13. 



Original design : Longerons — type 1 

Battens and diagonals — type 2 

Modified design: First bay — Longerons — type 1 

Battens and diagonals — type 2 

Other bays — Longerons — type 3 

Battens and diagonals — type 4 


Figure 13 
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STATIC LOADING 


The first problem cons 
load in the z direction at 
mary of the results. The t 
(point a) of the original s 

to design modifications is 

the corresponding w and 4> ^ 
than those for the original 
procedure with four global 
of the modified structure t 


idered is that of the static response due to a transverse 
the free end of the cantilever. Figure 14 shows a sum- 
ransverse displacement w and rotation at the free end 
tructure are given. The sensitivity of these quantities 

3X 8 

provided by {— } and {— «} . Also shown in Fig. 14 are 
° 9X 2 o 

for the modified structure (which are considerably larger 

structure) . The solution obtained using the proposed 

approximation vectors was identical to the direct solution 
o at least three significant digits. 


Static analysis 


Original 
structure 
X * 0 


Modified 
structure 
X - 1 


w at a 


<J >2 at a 


M 

0.102 

-0.0397 

9 u 

0.0737 

-0.0257 

) 

0.110 

-0.0388 




r - 2 

0. 394 

-0. 142 

r ■ 4 

0.394 

-0. 143 

Full system 

0.394 

-0.143 



Figure 14 






FREE VIBRATIONS 


The second problem considered is that of the free vibrations of the same lattice 
structure. Figure 15 shows a summary of the results. The first three eigenvalues 
(squares of the vibration frequencies) and their first two derivatives with respect 
to X are listed. The corresponding eigenvalues of the modified structure are also 
listed. The frequencies predicted by the proposed procedure with eight approximation 
vectors (four eigenvectors and their first derivatives with respect to X) and twelve 
approximation vectors (four eigenvectors and their first two derivatives with respect 
to X) are listed. The predictions of the eight-vector approximation are accurate for 
the first two eigenvalues, but not the succeeding ones. On the other hand, the pre- 
dictions of the twelve-vector approximation are accurate for the first three 
eigenvalues . 



Mode 

1 

2 

3 

Original 

10“ 

5 

P x Q 

0.324 

0.502 

1.815 

structure 


5 dfi 




X = 0 

10 

X dX 

-0.134 

-0. 207 

-0. 486 


10' 

5 & 2 n 
* ? 

-0.0911 

-0. 141 

-0.472 







Modified 

m‘5 v 

r = 8 

0.122 

0. 191 

0.383 

structure 

1U X 





X = 1 

Q 

r = 12 

0. 122 

0.191 

0.896 



Full 

0. 122 

0.191 

0.896 



system 





Figure 15 




SUMMARY 


In summary, a computational procedure has been developed for calculating the 
sensitivity derivatives of large structural systems as part of structural reanalysis 
(see Fig. 17). The three key elements of the procedure are: 

a) lumping of the large number of design variables into one (or small number 
of) tracing parameter (s) ; 

b) application of operator splitting/reduction technique; and, 

c) for very large problems use of multilevel substructuring technique. 

The proposed procedure can be considered as a general computational strategy for 
generating the response of the modified structure using large perturbations from the 
response of the original structure. 

For static problems the similarities between the proposed procedure and precon- 
ditioned conjugate gradient technique are identified and are exploited to provide a 
rational procedure for selecting the preconditioning matrix and a physical meaning 
for the preconditioned residual vectors. 

Future work includes : 

o extension to more complex structures and to shape design modifications 

o generation of sensitivity information with respect to design variables. 


• Computational procedure presented for calculating 
sensitivity derivatives as part of performing 
structural reanalysis for large-scale problems 

• Lumping of design variables into tracing parameter* s) 

• Application of operator splitting/reduction technique 

• Use of multilevel substructuring 

• Future work includes: 

• Extension to more complex structures and to shape 
design modifications 

• Generation of sensitivity information w.r.t. design variables 


Figure 17 
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